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Abstract 



The equilibrium properties of the Blume-Emery-Griffiths model with bilinear 
quenched disorder are studied for the case of attractive as well as repulsive 
biquadratic interactions. The global phase diagram of the system is calculated 
in the context of the replica symmetric mean field approximation. 



I. INTRODUCTION 

The Blume-Emery-Griffiths (BEG) [[L|] model has been successfully applied to explain the 
behavior of different physical systems such as He 3 -He 4 mixtures, microemulsions, semicon- 
ductor alloys, to quote only a few. In this paper we investigate the effect of quenched disorder 
on the mean field phase diagram of a version of the BEG model in which orientational and 
particle degrees of freedom are explicitly introduced. We consider the Hamiltonian: 

K 

' H X J 'i S ' S .i"'"j ~ I'Yl" ■ ' (!) 

i<j i<j i 

where Si = ±1, n, = 0, 1; the bilinear couplings J^- are quenched Gaussian random variables 
with average Jo/N and variance J 2 /N, and the sign of biquadratic interaction may be neg- 
ative. This Hamiltonian represents a general framework to study the complex behavior of 
many physical systems. In the simple case of pure systems it describes, for example, a fluid 
model with magnetic properties like those of polar liquids, which exhibits multiple exotic 
phase diagrams In presence of quenched disorder the Hamiltonian ([]]) implements 

different models whose rich phase diagrams are not yet fully explored. Some limiting cases 
include the standard Ising spin glass J7J (/i — > oo; rii = 1 Vz) and the Ghatak-Sherrington 
model |8|JT^J^- pT| , p!3| (K = 0). In the particular case of K = —1 we recover the Frustrated 



Ising Lattice Gas |H| , |ll)| that is a version of the Frustrated Percolation model [IB], whose 



properties suggest a possible close connection with the theories of structural glasses [fL7[] . Es- 



sentially, the model considers a lattice gas in a frustrated medium where the particles have 
an internal degree of freedom (given by its spin) that accounts, for example, for possible 
orientations of complex molecules in glass forming liquids. These steric effects are greatly 
responsible for the geometric frustration appearing in glass forming systems at low temper- 
atures or high densities. Besides that, the particles interact through a potential that may 
be attractive or repulsive depending on the value of K. 

Here we try to further elucidate the phase diagram properties of Hamiltonian ([!]) within 
a replica mean field theory approach. Specifically, the presence of disorder in the bilinear 
term leads to the appearance of a transition from a paramagnetic to a spin glass phase at 
low temperature or high density. When the biquadratic interaction is attractive or weakly 
repulsive we found, accordingly with the value of density, two spin glass transitions of differ- 
ent nature separated by a tricritical line: at high density the transition is continuous while at 
sufficiently low density the transition becomes discontinuous. When the particle repulsion is 
rather strong two new phases with a sub-lattice structure may appear: an antiquadrupolar 
phase and at lower temperatures, an antiquadrupolar spin glass phase. The antiquadrupolar 
spin glass transition line is always continuous whatever is the value of the chemical potential. 
When the bond weights are not symmetrically distributed, a ferromagnetic phase (as well as 
a ferrimagnetic one, depending on the value of Jo) may appear in the phase diagram. One 
can also identify a line where a dynamical instability appears, as in spin glass models with 



discontinuous transition like the Potts glass p8|,p!9| and the p-spin model 1 17,20]. Moreover, 



preliminary results on the stability analysis indicates the presence of a transition line from 
a 'phase' with one to a 'phase' with infinity replica symmetry breaking. This rich scenario 
and the unification picture provided by the model are the main motivations of this paper. 



II. THE PHASE DIAGRAM 

Many general properties of the phase diagram of Hamiltonian ([!]) may be observed by 
simply studying the replica symmetric solutions of the mean-field equations. In the following 
paragraphs we discuss how the many thermodynamic phases appear when the parameters 
of Hamiltonian (|l|) are changed. 

For the sake of clarity we work out separately the cases of attractive, K > 0, and 
repulsive, K < 0, biquadratic interaction. 

A. Attractive biquadratic interaction 

Using the replica method the free energy can be computed as: 

Pf = - lim 1 1 \n[Z n ] av , (2) 

n-»0 n 

where n is the numbers of replicas and [■ • -] av denotes the average over the disorder. We 
obtain: 



2 



/3/=lim- 

n-»0 n 



Z a<b Z a 



+ -(P 2 J 2 + 2[3K)J2d 2 a -\nZ' 



(3) 



where Z' is the single site partition function: 



Z' = Tr exp \ (3 2 J 2 V q ab S a n a S b n b + (3 V J m a 5 a n a 



rr 



(4) 



and the order parameters: density (d a ), magnetization (m a ), and Parisi overlap (q a b), are 
defined as 



d a = (n a ) , 
m a = (S a n a ) , 
q ab = (S a n a S b n l 



(5) 
(6) 
(7) 



In the replica symmetric approximation, m a = m, d a = d and q a b = <?(1 — Sab), one finds: 
f3f = l -p.J,m 2 - \f3 2 J 2 q 2 + \{f3 2 J 2 + 2/^)d 2 



Vz In 2 



1 + e" cosh(/3Jz v /g + (3 Join) 



where Vz = dze z2//2 is the standard Gaussian measure and we define: 

8 2 J 2 

E = ^-—(d-q)+(3Kd + (3n . 

The order parameters satisfy the saddle point equations: 

cosh(/3 J z^/q + [3 Join) 



d= Vz 



e~ s + cosh(/5Jz v /g + (3 Join) 

f smh. 2 ([3Jzyfq + (3J m) 

q = Vz 

m = J Vz 



[e H + cosh(/5 J z^Jq + (3 Join)} 2 
smh((3Jzy/q + (3 Join) 



e = + cosh(/3Jz v /g + (3 Join) 



(8) 
(9) 

(10) 

(11) 
(12) 



Notice that with a suitable transformation of our parameters (f3u — > (3u — In 2) we recover 
the equations obtained using spin-1 variables. 
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1. The case Jq = 



The simplest case in this situation corresponds to a zero disordered average, Jq = 0. As 
can be seen in Fig. [1], where the phase diagrams in the plane (T, /i) are depicted for several 
K, at low temperatures or high values of the chemical potential (high densities), there is a 
phase transition from a paramagnetic (P; m = q = 0) to a spin glass (SG; m = 0, q ^ 0) 
phase. For positive values of K (and Jo = 0), these are the only phases allowed to the 
system. They are separated by a continuous transition up to a tricritical point, T trc , below 
which a discontinuous transition is observed. We just note that the Sherrington-Kirkpatrick 
(SK) model with its critical temperature T c j J = 1 is recovered in two limits, both giving 
the highest possible density (d — 1): \i — > oo or K — > oo. 

We locate the continuous transition by expanding the saddle point equations for small 
values of q, what leads to 



J ( K J fi 

— — 1 + exp — - — — — — — 



T, 



J 2T r , T, 



(13) 



This equation is valid up to the tricritical point that may be located by expanding the single 
site partition function around the transition line where q a b = and d a = d f2lf . This leads 
to: 



- trc 

T 



T trc -3 + 2K/J+J9 + A(K/J) 2 -AK/J 



AK/J 



(14) 



17)/4 ~ 0.219 and fj, trc /J ~ -0.559 |4 



For instance, if K/ J = — 1 we get T trc / J — (5 - 
The line T trc (K) is also shown in Fig. |l|. 

The exploration of our equations in the T = limit offers a simple tool to understand 
the subtle properties of the phase diagrams reported above. Moreover, remembering that 
our model shares features with the frustrated percolation where no frustrated loops may be 
closed, it is interesting to study the behavior of the density d at T — 0. This behavior is 
depicted in Fig. pi. At T = the saddle point equations become: 



d = erfc 



where C = [3J{d — q) and: 



1 /l K 7 [I 
2dK2 C + 7 d+ l 



(15) 



C 




1 /I K , (JL 

-Yd[-2 C + 7 d+ l 



Actually, d satisfies these equations up to /x* given by 

J 



-K- 



(16) 



(17) 



below which the close packing configuration is never achieved (d < 1); above //*, d = 1 and 
C = y2/n. The point //* (see Fig. ||) individuates a clear cusp in the behavior of d(fi)\ T=0 
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and seems to be a characteristic value for the system. We also observe that for high values 



of K, K > y2/7r, the region with < d < 1 disappears. 

From the tricritical point other two lines depart in the plane (/i,T) (see Fig.[TJ). They 
correspond to a thermodynamic line of discontinuous transitions from a paramagnetic to a 
SG phase, located where the free energy of the paramagnetic and spin glass solution are 
equal; and a purely dynamical transition line whose presence is due to the large number of 
TAP metastable states |22|] which can trap the system for very long times (infinite time in 
the mean field model) The study of this metastable glassy states within a suitable 



generalization of the TAP approach to the present model is in progress [23 



In order to identify the dynamical transition line we should compute the smallest eigen- 



value (replicon) of the 1RSB stability matrix and impose the marginality condition |26 



An alternative procedure amounts to expand the 1RSB free energy (see appendix) around 
the point m = 1, where m denotes here the size of diagonal blocks in the Parisi ansatz 
T9| . However at level of the approximation we work (replica symmetry) we can guess that 



the dynamical transition line is located near the line where a non zero solution of replica 
symmetric saddle point equations first appears; however, to check this point a more detailed 
investigation will be need. 

When the temperature and chemical potential are close enough to the discontinuous 
transition line we expect that the glassy phase of the system is exactly described by only 
one step of replica symmetry breaking [p3fl . On the other hand, when — > oo (the SK 
model), the full RSB scheme, with infinite steps, is needed in order to obtain the correct 
answer. Thus, a further 'replica transition line' from a SG phase with one (or finite) replica 
symmetry breaking level to a SG phase with infinite levels of replica breaking is expected to 
start at the tricritical point and terminate on the T = axis in some characteristic point, 
probably near fi*. This transition line is expected to separate qualitatively different aging 
behaviors of the system [27|-p9|. 

We can also apply a magnetic field to the system by including a term like —hJ^iSifii 
to the Hamiltonian eq. QT|) (this is done by adding the term (3m to the argument of the 
hyperbolic functions in the saddle point equations). The continuous transition is destroyed 
as soon as h ^ and the first order transition ends in a critical point, the system presenting 
the characteristic wing shaped phase diagram. The effects of a magnetic field, as well as 
the stability of the RS solution in its presence were extensively studied for the Ghatak- 
Sherrington model ||13|. 



2. The case J ^ 

Having given some details for the case Jo — 0, we face below the more general situation 
when J is different from zero ||32|| . There are now three kind of phases (see Fig. [3D: a para- 



magnetic and a spin glass phase with essentially the same characteristics we have discussed 
above, along with a ferromagnetic phase (F; g ^ 0, m ^ 0) which appears, as usual, in the 
low temperature region (or high fi) when Jo ^> J. As in the SK model, this ferromagnetic 
phase is reentrant (in RS), although here the degree of reentrance may vary. 

When the transitions between these phases are continuous (e.g., Fig.^J) , the boundary 
lines may be analytically obtained. The boundary between F and P is given by 
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while the boundary between SG and F obeys 



T c = J (d - q) , 



(19) 



where d and q are obtained solving the saddle point equations with m — 0. And the boundary 
between SG and P is our former result eq. flT3"|) . Notice that in the limit [i (or K)—> oo , 
we recover the SK boundaries (T c = J for SG-P, T c = J for F-P and T c = J (l — g) for 
SG-F). 

Changing the value of the chemical potential the transitions may become discontinuous, 
as can be seen in Fig. |]for \i = —0.6 and K = —1, where a tricritical point and an endpoint 
show up. 

As usual the analysis at T = is interesting. In this limit we get the following equations: 



d 



m 



-erfc 



■hC — Kd — \i — J^ra 



2d 



1 f (-\C-Kd- n + hm 
H — erfc — = 

2 V V2d 



(20) 



-erfc 



\C — Kd — ji — Jom" 



2d 



1 / — \C — Kd — ji + J m" 
-erfc — = 

2 V V2d 



(21) 



and 



C 




cosh 



(lC + Kd + fi) +Jlm 2 
2d 



J ° m ( C , ISA , 

— [2 +Kd+ J 



(22) 



We can then individuate the location of the continuous transition line, in the plane 
(Jo,/i), between the F and SG phases. It occurs where m — > so, expanding the above 
equation in m, we get: 



Jo 




cxp 



[\C + Kd + fi 
2d 



(23) 



Here d and C are evaluated in the points where m = 0. This equation is valid up to a given 
jj,* = —K — 1/ \[2~tx (note that since m = 0on this transition line, we just recover our former 

result) and above this value, we find J = yV/2 ~ 1.253. 

As is well known in the SK model, when J grows too much, the SG phase disappears. 
A similar phenomenon happens here, but the order of the transition between the phases F 
and P depends on the value of fi, as expected. Notice that we are considering K — — 1, 



6 



but in the next section we will show that for this value there are no antiquadrupolar (glassy 
or not) phases and the characteristics presented here are also general for greater values of 
K. We expect that for negative values of K or J , there is also the possibility of having 
ferromagnetic or anti-ferromagnetic ordering, respectively The former is characterized by 
having different magnetizations in the sublattices (itla ^ itlb) and the latter by having 
opposite magnetizations (m^ = -ms). The analysis of these phases is beyond the scope of 
this paper. 

B. Repulsive biquadratic interaction 

Because the sign of the biquadratic coupling is now negative we have to take into account 
the possibility of further phase ordering within a two sub-lattice structure. This for instance 
is the case studied in refs. |2|-|6| in absence of quenched disorder and in ref. but only 
for fi = 0. When oo we recover the anti- ferromagnetic Sherrington-Kirkpatrick model, 
studied in ref. |34| , and following their prescription, the Hamiltonian can be rewritten as: 



E E< - 



* 3 



i -J 



(24) 



a=A,B i 



where a = A, B is the index of the two sub-lattices. Using the replica method we obtain 
the following free energy: 



/?/ = lim- 

n—>Q n 



\p 2 J 2 E qaQb + h Jo E ™>£ 

/ a<b A a 



+\((3 2 J 2 + 2(3K) ^ d\d B - i In Z A Z B 



(25) 



where Z a (a = A, B) is the single site partition function of the sub-lattice a: 



Z a = Tr exp < (3 2 J 2 V q ab S a n a S b n b + (3 V J m a S a 

S ' U I a<b 



n 



+E 



'(3 2 J 2 



+ pK)d a + pfi 



(26) 



In the replica symmetric approximation the free energy can be easily computed: 
Pf = lf3J m A m B - ^(3 2 J 2 q A q B + ^((3 2 J 2 + 2(3K)d A d B 
-]- J2 f^zln2 1 + e Ea cosh({3 J Zy/q^ + f3J m c 

a=A,B 



where 



f3 2 J 2 



{d a - q a ) + f3Kd a + (3 fx . 



(27) 
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From the free energy, one derives the replica symmetric saddle point equations: 



Vz cosh((3 J Zy/q^ + f3J m B ) 

e~ SB + cosh(f3Jzy/qB + f3J rriB) ' 

U ~ 1 Z [e- s B + CO sh(f3Jz^ + (3J m B )} 2 ' 1 ' 



smh(/5Jz v ^i + 0J o m B ) 

m A = / Vz — = tttt? — 1= — 7Tr r • (30) 

J e~- B + cosh(/? JZy/q^ + (jJ m B ) 

The analogous equations for d B , q B and m B are obtained from the previous one exchanging 
A <-> £>. A representative case is shown in Fig. || where both gA(s) and g?a(B) are presented 
in the SG and AG phases (see definitions below). 

As stated above, new different phases emerge here. At high temperatures (or low /x), 
when the orientational degrees of freedom are not interacting (q A — q B — 0), there are 
two possible orderings: a paramagnetic (or quadrupolar) phase, P ,with d A = d B ^ and 
an antiquadrupolar one, AQ, with c/^ ^ d B . The name quadrupolar is reminiscent from 
the spin-1 representation (r, = ^rij) where it labels the ordering of variables rf = rii. At 
low temperatures (or high /i), the orientational degrees of freedom become important and 
different glassy phases show up. First, there is a spin glass phase (SG) with d A = d B and 
<?a = qB 7^ 0). This is the only glassy phase if the effects of particle repulsion are not very 
strong (K not too negative). On the other side, when K is highly negative, besides the SG 
phase, there is also an antiquadrupolar glass phase, AG, where the sub-lattice symmetry is 
broken (d A ^ d B , q A ^ q B ). 

As in the previous section, when J is allowed to be nonzero, the system develops a mag- 
netization and additional, ordered phases may appear: when J ^> J at low temperatures 
(or high /i) a ferromagnetic phase (F; m A = m B ^ 0, d A = d B , q A = q B ^ 0) is encountered 
if K is not too negative, or either a ferrimagnetic phase is entered (J; m A ^ m B , d A ^ d B , 
qA 7^ qB)- Also, with J —J, an anti-ferromagnetic phase {m A = —m B ^ 0) may shows 
up. 

From now on we consider Jo = 0. As we decrease K from a positive value, there is 
a point, K AG ~ —1.46, where the AG phase first appears, growing inside the previously 
described SG phase (see Fig. [?[). At a characteristic value, K A q = —3/2 — y/2 ~ —2.91, 
an antiquadrupolar phase appears between the paramagnetic, the SG and the AG phases, 
as depicted in Fig. [5[ This value is obtained noticing that there are two points where the 
P-AQ line intercepts the AQ-SG line (see Fig. ||), given by 




T ± = ^? k +t;±\K 2 + 3K + - . (31) 



When both are equal, we get K A q. As soon as the AQ phase appears, the SG phase is 
divided in two regions. At even lower values of K (see Fig. the AQ grows, while the left 
SG phase shrinks and the right one moves to higher values of /i. This SG phase disappears 
as K — > — oo. We can also notice that the AG phase is invaded by the AQ as we approach 
this limit, and eventually we have only the P and AQ phases, as we could expect since in 
this limit we have a lattice gas with repulsive interactions. 



S 



With the same techniques adopted in the previous paragraph, we can analytically locate 
the P-AQ continuous transition line reported in Figs. || and |9|: 



T H d(d-l) 



where d satisfies the equation: 



dexp 



1 - d 



[1 — d) exp(/3^) . 



Analogously, the AQ-AG continuous transition line is given by: 

— — dj^ds , 



(32) 



(33) 



(34) 



where dA and ds satisfy the saddle point equations with qA = Qb = 0. 

Finally, the AG-SG transition line can be found expanding the saddle point equations 
in the small quantities qA — gU — d B , representing the staggered order parameter. 

The stability of the solutions presented here can be studied by evaluating the eigenvalues 
of the matrix of gaussian fluctuations. The eigenvalues structure is analogous to the one 



found in the SK model ||30| , |3T| with non zero magnetic field since here the chemical potential 
plays a similar role [ 14j, 32 ,[33 . However, in this case the eigenvalues may be negative or even 
complex 



, |T2| , |3"3"|j . The numerical study of their behavior |33| shows that the paramagnetic 
and antiquadrupolar phases are inside the stability region while both glass phases (AG 
and SG) are unstable. The critical line P-AQ, P-SG and AQ-AG are found to be in the 
stability region, their location depending on the chemical potential, and should not change 
when replica symmetry is broken; while the transition line AG-SG is found to be completely 
inside the instability region and its precise location, beyond the aims of this paper, may be 
found only within the Parisi scheme of RSB. 



III. CONCLUSIONS 

We have studied the global phase diagram of a spin glass version of the Blume-Emery- 
Grifliths model where the bilinear couplings are quenched gaussian random variables. The 
particles, besides the steric effects due to complex molecular structure (represented by the 
possible spin orientations), are subject to a potential that may be either attractive or repul- 
sive. These steric effects that prevent close packing configuration are common in materials 
like glasses and poured sand. 

This model displays a large variety of interesting critical behaviors. When the particles 
interaction is attractive or weakly repulsive the transition between the paramagnetic and 
spin glass phase may be either continuous or discontinuous, depending on the value of the 
chemical potential; these different behaviors are separated by a tricritical line. For strong 
repulsive interaction between particles, new different phases with a sub-lattice structure 
emerge: the antiquadrupolar and the antiquadrupolar glassy phase. These phases seems to 
be separated by a continuous transition whatever is the value of the chemical potential (at 
least in the range of K studied here). When the quenched disorder is not symmetrically 
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distributed, also a ferromagnetic phase or a ferrimagnetic phase may appear in the phase 
diagram. 

A further transition line where a dynamical instability appears may also be found as 
happens in other spin glass models with discontinuous transition. Indeed, many interest- 
ing properties emerge if we consider the dynamical behavior of systems whose equilibrium 
properties are described by Hamiltonian (JX|) . For example, the model with diffusive particle 
dynamics shows strong glassy behavior characterized by diverging relaxation times, vanish- 
ing diffusivity and breakdown of the Debye-Stokes-Einstein law [O]; while the spin-1 version 



of the model with Monte Carlo dynamics displays different aging regimes due to the existence 
of spin-glass transitions of different nature. In the region where the spin-glass transition is 
discontinuous the aging behavior is stable against non-relaxational perturbation |35j and 
quite similar to the one observed in non-relaxational dynamics of the spherical p-spin glass 
Furthermore when a gravitational term is added to the Hamiltonian ([I]), strong 



model 36 



links appears with granular media and their complex dynamical behavior as logarithmic 



compaction [37 
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APPENDIX A: THE 1RSB SOLUTION 

In this appendix we present the solution of the Blume-Emery-Grifhths spin-glass model 
with one-step of replica symmetry breaking (1RSB). The following relations are the starting 
point to get the phase diagram corrections and the dynamical transition line. Without loss 
of generality we consider the case Jo = where the magnetization is zero. Following the 
Parisi scheme, the n replicas are divided in n/m blocks containing m replicas. Different 
replicas in the same block have overlap qi while those in different blocks have overlap go- 
Thus, for the case of attractive particles interaction (K > 0), the 1RSB free energy reads: 



1 



f3f x = --^J 2 [(1 - m)q( + mql - tf J + -fiKd* - In 2 
1 



1 



Vz In / Vz x (1 + e =1 cosh Qi) m (Al) 



where: 



fii = (3Jz ^ + z l ^Jq l - q (A2) 
2i = ^~(d - qi) + P(jjl + Kd) (A3) 



The saddle point equations are: 
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, ^ / e =1 cosh Qi \ . . , . 

d =l^ TT^ihk (A4) 



-2 



* = l v * i + eH.coshnJ < A5) 



/"^ / e^sinhf^ \ . 
■"=/^( 1 + e H, cosh 1 ti J (A6) 

(d> qi > go) an d m satisfies the equation: 

^m 2 p 2 J 2 (ql - gg) - m ^ hi(l + e Hl coshfti) 

+ J Vz Q In J V Zl (1 + e Hl cosh fix)" 1 = (A7) 

Here the overbar denotes the average: 

= JV Zl (l + e^ cosh^rX 

/Psi (1 + e s icoshfi 1 ) m 1 ' 

Following the same procedure, for the case of repulsive particles interaction (K < 0) we 
obtain: 

Pfi = - \(3 2j2 K 1 - m)q lA q lB + mq 0A q 0B ] (A9) 
-- A [f3 2 J 2 + 2pK]d 2 -\n2 

~ ^ E fl / In / Z>Z! A™(z , (A10) 



where: 



A*(^o,2i) = l + e =lQ cosh(/5J2;o v /go^ + (3Jzi^fqi a - q 0a ) 
P 2 J 2 



— la 



■(d-9ia) +P(pL + Kd) 



2 

It's straightforward, but tedious exercise, to obtain the saddle point equations. 
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FIGURES 



1.0 




FIG. 1. Phase diagram for several values of K/J showing the paramagnetic (above) and spin 
glass phase (below the line). The continuous line stands for the second order transition while 
the long dash line is the discontinuous transition. Also depicted (dotted line) the line where the 
solution q ^ first appears; while the small dash line is a line of tricritical points. 
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FIG. 2. The order parameter d at T = for K/J = — 1,0 and 1. The dotted part signals the 
region where non zero solution of saddle point equations first appears. Observe that there is a 
particular value of the chemical potential, /j,*, where the density becomes lower than 1. 
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Jo 

FIG. 3. Phase diagram for K = — 1 and /U = 0. The reentrance is an artifact of the RS 
approximation. All transition lines are continuous and look much the same as the SK boundaries. 
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FIG. 4. Phase diagram for K = — 1 and ji = —0.6. The continuous lines are continuous 
transitions while the dashed lines stand for discontinuous transitions. Notice that the line SG-F 
ends in an endpoint and the transition F-P is discontinuous up to the tricritical point. 
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FIG. 5. Phase diagram at T = for K 
meet. 



-1 showing the critical endpoint when both curves 
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FIG. 6. Plot of the order parameters q (solid line) and d (dashed line) as a function of the 
chemical potential \i for K/J = —3 and a fixed temperature (T = 0.1). Both SG and AG phases 
are presented. In particular, the AG solution is the region were both q and d are twofolded. This 
figure is an horizontal cross section of the fig. ||. 
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